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Abstract. We report on a comparative study of the valence band electronic 
structure of NiO as bulk material and of NiO as impurity in MgO. From the 
impurity we have been able to determine reliably the parameters which describe 
the local correlations, thereby establishing the compensated-spin character of the 
first ionization state or the state created by hole doping. Using bulk-sensitive 
x-ray photoemission we identify pronounced satellite features in the valence band 
of bulk NiO which cannot be explained by single-site many body approaches nor 
by mean field calculations. We infer the presence of screening processes involving 
local quasi-core states in the valence band and non-local coherent many body 
states. These processes are strong and the propagation of an extra hole in the 
valence band of NiO will therefore be accompanied by a range of high energy 
excitations. This in turn will make the observation of the dispersion relations in 
the Ni 3d bands difficult, also because the effective band width is no more than 
0.25 eV as estimated from multi-site calculations. 
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NiO is a benchmark system in solid state physics. It crystalhzes in the NaCl 
structure, has a partially filled 3d shell (Ni^+ d^), and is an antiferromagnetic insulator 
with a Neel temperature of 523 K [I]. It was pointed out early on by de Boer and 
Verwey [5] that many of the properties of the 3d transition metal compounds do not 
agree with the predictions of band theory, e.g., standard band theory predicts NiO to 
be metallic. A qualitative explanation was proposed in terms of the Mott-Hubbard 
model [21 S] in which the on-site Ni 3d-3c? Coulomb interaction plays a decisive role. 

An early ab initio attempt to fix the shortcoming of band theory was to treat NiO 
as a Slater insulator in which the doubling of the unit cell allows for the existence of 
a gap [3 mi?]. However, the calculated gap of about 0.2 eV [5] turned out to be much 
too small: A combined photoemission (PES) and bremsstrahlung-isochromat (BIS) 
spectroscopy study showed that the band gap is 4.3 eV [5] and established thereby 
the correlated nature of NiO. The inclusion of a self-interaction-correction (SIC) or 
Hubbard U term to the density-functional formalism may provide a justification for 
the magnitude of the experimental band gap [H [TU] . 

Yet, one of the most direct methods to critically test the accuracy of the different 
approaches, is to determine the excitation spectrum associated with the introduction 
of an extra particle into the system [13]. Curve (a) in Fig. 1 displays the valence 
band x-ray photoemission spectrum (XPS, hv = 1486.6 eV) of an in situ cleaved 
NiO single crystal. This spectrum represents essentially the Ni 3d spectral weight 
since the photoionization cross section of the O 2p is relatively small [M]. One can 
clearly observe from curve (e) in Fig. 1 that the Ni 3d density of states calculated 
by band theory (in the local density approximation, LDA) does not match at all: It 
has a Fermi cut-off and the line shape is completely different. The inclusion of the 
Hubbard U in the calculations (LDA-f U) does not solve the line shape problem, see 
curve (d). All this demonstrates the shortcomings of mean field theories to describe 
spectra associated with the fundamental one-particle Green's function of the system 

[iniiii]. 

A completely different approach is to give up the translational symmetry of the 
system in order to focus on the local correlations and, especially, the dynamics of the 
propagation of the injected particle. Curve (c) of Fig. 1 shows the Ni 3d spectral weight 
from an early cluster configuration-interaction calculation by Fujimori and Minami 
[12], which also includes the full atomic multiplet theory. The agreement with the 
experimental spectrum is extremely good. Nevertheless, a later cluster calculation by 
van Elp et al. [11] arrived at a less satisfactory result: Peak B has almost disappeared 
in the calculation, see curve (b). The prime motivation to use a different set of 
model parameters here is to infer that the first ionization state is low spin i^E) [16] 
rather than the Hund's rule high spin (^T), analo gous to the case of Zhang-Rice 
singlets in the cuprates [ITl [18] . Recent developments combining LDA with dynamical 
mean field [TU |20l |2T1 [22] or GW approaches [23] yield Ni 3d spectral weights which 
deviate in important details from the experimental spectrum. These discrepancies 
between the experiment and the later theoretical simulations [TTl [191 HOI [HI [22] do 
not provide confidence that one has made progress in understanding the nature of the 
first ionization state. 

The issues that we need to address now are threefold. First of all we have to 
establish whether the XPS valence band spectrum in Fig. 1 is truly representative 
for bulk NiO. There are reports in the literature claiming that certain satellite peaks 
in the Ni 2p spectrum are due to surface effects [M] |25l [26l [27] . Second, we have 
to determine to what extent a single-site many body approach can be utilized to 
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Figure 1. (Color online) Valence band XPS (1486.6 eV) spectrum of an in 
situ cleaved NiO single crystal. The results of two single-site cluster calculations 
(reproduced from Refs. '11' and '12'), and LDA and LDA+U calculations are also 
included for comparison. 

describe the electronic structure of NiO for which band formation is also essential. 
Third, we need to identify the nature of the first ionization state in the framework of 
a local ansatz. To this end we measured the valence band of NiO utilizing the more 
bulk-sensitive hard x-ray photoelectron spectroscopy (HAXPES) and we investigated 
experimentally the electronic structure of NiO impurities in MgO. 

The XPS data {hv = 1486.6 eV) on in situ cleaved NiO single crystals were 
recorded using a Vacuum Generators twin crystal monochromator A\-Ka source and 
an Scienta SES-100 analyzer, with an overall energy resolution set to 0.35 eV. The 
HAXPES data {hiy = 6500 eV) were taken at the Taiwan beamline BL12XU of SPring- 
8 in Hyogo, Japan using an MBS A-IHE analyzer. The overall energy resolution was 
set to 0.35 eV. The NiO impurity in MgO system was prepared in situ as 10-20 nm thin 
films on polycrystalline Ag by means of molecular beam epitaxy. The measurements 
were performed at the llAl Dragon beamline of the NSRRC in Hsinchu, Taiwan. The 
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Figure 2. (Color online) Valence band photoemission spectra of an in situ cleaved 
NiO single crystal recorded using 1486.6 eV (XPS) and 6500 eV (HAXPES) 
photons. 



photoemission spectra were recorded at the Cooper minimum of Ag Ad {hv = 140 eV) 
[28j using a Scienta SES-100 analyzer with an overah energy resolution set at about 
0.15 eV. 

In Fig. [5] we show the valence band photoemission spectra of a freshly cleaved NiO 
bulk crystal, taken with a photon energy of 1486.6 eV (XPS) and 6500 eV (HAXPES). 
By increasing the photon energy we increase also the kinetic energy of the outgoing 
photoelectron and, thus, also the inelastic mean free path. One can estimate that the 
probing depth is then enhanced from about 15 A to roughly 80 A [2^. We observe 
that the spectra are very similar. We, thus, conclude that the XPS data as displayed 
in Figs. 1 and 2 is representative for the NiO bulk material and that the contribution 
of surface effects [24j [25l [26l [27] can be safely neglected. To be specific: Peak B is 
intrinsic for bulk NiO. We would like to note that increasing the photon energy from 
1486.6 eV to 6500 eV does not alter much the Ni Zd character of the spectrum. The 
O 2p photoionization cross section relative to that of the Ni Sd remains very small, it 
changes from 1/13 to only 1/10 [13], meaning that peak B truly belongs to the Ni id 
spectral weight and not to the O 2p [21]. 

The valence band spectrum of the Nio.05Mgo.95O impurity system is shown in 
Fig. |3] together with the spectrum of an MgO reference thin film grown simultaneously 
under identical oxygen and substrate conditions. The Nio.05Mgo.95O film (and also 
the MgO film) was capped by 2 monolayers of MgO in order to prevent the surface 
termination to have an effect on the local electronic structure of the Ni impurity. The 
inset in the figure displays the Ni ^2,3 x-ray absorption spectra of the Nio.05Mgo.95O 
and the NiO bulk. The spectra are essentially identical, verifying that the Ni in the 
Nio.05Mgo.95O has very similar local surrounding (NiOe octahedra) as in the bulk. 

The valence band spectra of the Nio.05Mgo.95O and MgO systems are normalized 
to their O 2s core level intensities. Both are dominated by the O 2p valence band, yet, 
there are clear differences between them due to the presence or absence of the 5% NiO 
impurity. The difference spectrum multiplied by a factor of 6 is given by the red curve 
in Fig. [3j The line shape remains the same for films with lower Ni concentrations, 
as is the case for that of the Ni 2p [3^. This curve represents essentially the Ni Sd 
spectral weight of the NiO impurity since the photoionization cross section of the Ni 
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Figure 3. (Color online) Extraction of the NiO impurity valence band 
photoemission spectrum: Valence band spectra of Nio.05Mgo.95O and an MgO 
reference, together with the resulting difference spectrum. Also included is the 
result of a single-site NiOg configuration-interaction cluster calculation. The inset 
shows the Ni L2,3 x-ray absorption spectrum of the Nio.05Mgo.95O together with 
that of bulk NiO. 



3d is an order of magnitude larger than that of the O 2p at the photon energy used 
[M]. Remarkable is that it is different from the spectrum of bulk NiO as shown in 
Figs. 1 and 2. The impurity spectrum lacks specifically peak B which is prominently 
present in the bulk spectrum. 

To interpret and understand the impurity spectrum, we have performed 
simulations using the well-proven configuration-interaction cluster model which 
includes the full atomic multiplet theory [211 [321 [33] ■ The simulations have been 
carried out for a NiOg cluster using the program XTLS 8.3 [3T] . 

The bottom curve in Fig. [3] shows the Ni 3d one-electron removal spectrum from 
the cluster calculation. The agreement with the experiment is very satisfactory. In 
order to achieve this, we have started the calculations by using parameter values which 
were suggested from earlier studies on NiO [TTl [311 [311 [33 ■ We then fine-tune the 
parameters describing the octahedral crystal and ligand fields, and also the difference 
between the Hubbard U and the O 2p-Ni 3c? charge transfer energy [36]. The crucial 
issue here is to obtain a main line (peak A) without having another feature appearing 
at about 2 eV higher energies (peak B) as was the case in the simulations by Fujimori 
and Minami [12|] and by van Elp et al. [Tl]. This has implications for the energetics 
of the states making up the valence band as we explain in the following. 

A detailed look at the cluster calculations displayed in Fig. 1 shows that peak A 
is given by the "'Ti final state of the Ni Sd^ multiplet structure while peak B is due to 
the ^Ti. Avoiding the appearance of peak B means that the energy splitting between 
these two states must be made smaller, e.g., 1 eV or less. This is what we have done 
in our simulation in Fig. 3, using different but equally reasonable parameter values 
[36] . The consequences for the physics are, yet, quite far reaching. Given the fact that 
various x-ray absorption studies find an effective octahedral crystal and ligand field 
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Figure 4. (Color online) Comparison of the valence band photoemission spectra 
of bulk NiO and NiO impurity in MgO. Also included is the simulated Ni 3d and 
O 2p spectral weights of the NiO valence band from a LDA+DMFT calculation 
(reproduced from Ref. 1211 ) and the result of a three-site NiaOie configuration- 
interaction cluster calculation. 

splitting of about 1.65 eV [311[3S], i.e., the splitting between the isospin ^Ti and 
states, we arrive at the conclusion that the must be lower in energy than the *Ti 
by 0.65 eV or more. This is what we read from our results in Fig. 3. In other words, 
our impurity study provides the spectroscopic evidence that the first ionization state 
has a compensated-spin character rather than the Hund's rule high-spin. This in turn 
justifies that the ground state of a hole doped NiO system may indeed be low-spin in 
nature [IB] . 

We now return to the problem of the bulk NiO valence band spectrum. Fig. |4] 
shows the Ni 3d spectral weight taken with XPS and compares it with the spectra of 
the NiO impurity and of the single-site LDA+DMFT calculation [21] . One can clearly 
observe that peak B is absent in the impurity as well as in the Ni 3d spectral weight of 
the single-site calculation. In fact, one could infer that the calculation reproduces quite 
well the impurity spectrum, with perhaps some discrepancies due to the incomplete 
implementation of the multiplct structure of the on-site Coulomb interactions. Yet, 
the discrepancy with the bulk spectrum strongly suggests that the origin of peak B 
must be sought in non-local correlations, i.e., effects which cannot be included in a 
single-site approach. 

Our suggestion is that peak B is due to non-local screening processes involving the 
formation of low-energetic coherent many body states on neighboring NiO clusters, 
which are of the ^E type as we have shown above. The mechanism is analogous as 
proposed earlier for the Ni 2p core level spectrum of bulk NiO [37] , but the application 
of it for the valence band is only valid for local states which are relatively stable against 
band formation. This may not be applicable for the '^E state, which is a state in which 
a hole is injected in an eg orbital starting from the 3d^ ground state [35]. This 
hole can be expected to readily propagate in the lattice since the hopping between the 
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Figure 5. (Color online) Non-local screening in valence band photoemission 
on NiO: (a) Creation of the atomic-like (quasi-core) ^Ti hole state by the 
photoemission process, (b) Screening by a next nearest neighbor NiOe cluster 
producing a coherent low-energetic hole state there. 

Ni 3d(eg) and O 2p{a) orbitals are rather large [TTl |37], yet, it may leave behind an 
energetically costly wake of wrong spins in the antiferromagnetic lattice. In any case, 
it would not be meaningful to describe its band formation as a low-energy screening 
process involving neighboring "^E states fSH^. 

However, for the main peak of the bulk NiO spectrum, i.e., the ^Ti state, we infer 
that we can make a meaningful approximation by using the coherent screening 
model. The ''Ti consists of a hole injected into the t2g orbital, and its ability to move 
is rather limited since the overlap between the Ni 3(i(i2g) and O 2p(7r) is small. One 
could consider the ^Ti as a localized quasi-core state. We then can invoke the non-local 
screening process as follows: After the creation of the ^Ti state, an Cg electron from 
a neighboring NiO cluster hops onto the Ni site, leaving behind a coherent ^E hole 
state on that neighbor. A sketch for this process is given in Fig. [5] These two states 
are energetically almost degenerate [37|, and the Ni id{eg) and O 2p((T) hybridization 
between them is then strong enough to produce two peaks: Not only the main peak 
A but also the satellite peak B. 

To confirm our assignments, we have performed a NiaOig cluster calculation 
consisting of three edge-shared NiOg octahedra. While all the O 2p and Ni Sd orbitals 
are included for the NiOg octahedron in the center where the photo-excitation takes 
place, those on the other parts of the cluster are replaced by a reduced basis set using 
the method in Ref. [39] . The results are displayed in Fig. |4] and demonstrate the 
presence of both peaks A and B. Note that we have used the same parameters as for 
the single-site calculation which produces only peak A, see Fig. [S] We also note that 
the energy difference between peaks A and B is somewhat smaller and the intensity 
of peak B is slightly larger than those of the experiments. This can be explained by 
the fact that the number of neighboring Ni sites is only two in the NisOig cluster: the 
energy difference will increase and the intensity of peak B will decrease for a larger 
number of neighboring sites [39] . 
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To summarize: We have succeeded to determine reliably the Ni Sd valence 
band spectra representative for bulk NiO as well as for NiO as an impurity system. 
From the impurity data we are able to extract the local electronic structure and the 
correlations herein, thereby establishing firmly the compensated-spin character of the 
first ionization state. By comparing the bulk with the impurity system, we were 
able to identify features in the bulk NiO spectrum which are caused by screening 
processes involving local quasi-core valence band states and non-local low-energetic 
many body states. These processes are strong and cause considerable redistribution 
of the spectral weight when an extra hole is injected into NiO. This explains why 
fine features related to the dispersion of the Ni 3d band is difficult to detect [40] , also 
since the effective band width is only of order 0.25 eV as estimated from our multi-site 
cluster calculations. 
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